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Abstract 

We present results of a Monte Carlo study of temperature-programmed des¬ 
orption in a model system with attractive lateral interactions. It is shown 
that even for weak interactions there are large shifts of the peak maximum 
temperatures with initial coverage. The system has a transition temperature 
below which the desorption has a negative order. An analytical expression for 
this temperature is derived. The relation between the model and real systems 
is discussed. 
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I. INTRODUCTION. 


Many studies have been published on Monte Carlo (MC) simulation of desorption from 
single-crystal surfaces, M but little attention has been given to the effect of attractive lat¬ 
eral interactions. There may be three reasons for this. First, temperature-programmed 
desorption (TPD) spectra with attractive lateral interactions appear less interesting than 
TPD spectra with repulsive lateral interactions, which can show multiple peaks even with 
just one adsorbate. Second, repulsive lateral interactions are more common than attractive 
lateral interactions. This can be understood using the bond-order conservation principled 
Nearest neighbor interactions are often repulsive, because nearest neighbors bond to the 
same substrate atom(s). Multiple bonds of a substrate atom are generally weaker than 
when a substrate atom bonds just to one adsorbate, so it’s energetically more favorable for 
adsorbates to stay apart. Only when the direct adsorbate-adsorbate interaction is strong 
with respect to the substrate-adsorbate interaction then one will find attractive lateral in¬ 
teractions. This is the case for physisorbed adsorbates. Third, a good quantitative MC 
method is necessary to understand the effects of attractive lateral interactions. In particu¬ 
lar, it is necessary to have an MC method with the correct real time, instead of MC time, 
dependence. Such methods have only recently been developed.!0 It is known that attrac¬ 
tive lateral interactions cause peaks in TPD spectra to shift to higher temperature and to 
become smaller .011 In this paper we will concentrate on shifts of peaks with initial coverage. 


II. COMPUTATIONAL DETAILS. 

We have used MC modeling of the reactions based on the master equation 

dP 

-J .- = £ \ w «e p i< - w 0° p «l (!) 

ai p 

where P Q is the probability of having the adlayer in configuration a. The IP’s are transition 
probabilities per unit time corresponding to various reactions. The master equation is solved 
using the first-reaction method in which a time for each reaction to occur is generated from 
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an appropriate probability distribution, and then the adlayer is modified according to the 
first reaction after which the procedure is repeated.0 Details of this method can be found 
elsewhere .0 The transition probabilities per unit time are given by 

E a g 


W a g = u a g exp 


knT 


( 2 ) 


where E a g and u a g are the activation energy and prefactor, respectively, of the reaction (3 —> 
a that transforms conhguration (3 into configuration a. The W s may be time-dependent 
through the temperature. For the simulation of TPD spectra 


T = T 0 + Bt 


( 3 ) 


holds, where T 0 is the temperature at the begin of the simulation, t is time, and B is the 
(constant) heating rate. 

For all the desorption reactions we have used the same prefactor u a g = uq. The activation 
energy is given by 


E a g — E 0 (l + 'jNag), (4) 

where E 0 is the activation energy for an isolated adsorbate, N a g is the number of (nearest) 
neighbors of the desorbing adsorbate in the reaction f3 —> a, and 7 is a positive dimensionless 
parameters determining the strength of the interactions between the adsorbates. The linear 
dependence of the activation energy on the number of neighbors seems to be appropriate 
for the systems that are best known for attractive adsorbate-adsorbate interaction, namely 
noble gas atoms, but we will also look at other dependencies on the number of neighbors 
later on. 

Apart from desorption we also have diffusion. This is modeled by jumps of adsorbates to 
neighboring sites. The diffusion should be fast with respect to desorption, because we want 
a thermally equilibrated adlayer before each desorption. We found that for values 7 < 0.1 a 
prefactor for diffusion of v a g = 20v 0 was sufficient. We used the same expression ([|) for the 
activation energy as for the desorption. The number of neighbors in that expression refers 
to the situation prior to the jump of the adsorbate. 
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Most MC simulations were done on a square lattice of dimensions 100 x 100. This 
effectively models an infinite lattice. Although we didn’t do extensive simulations of lattices 
with different numbers of neighboring sites Z, we will present an analysis that holds for any 
value of Z. All simulations have been done with our code PIZZAZZ.0 

III. RESULTS AND DISCUSSION. 

The MC simulations that have just been described depend on five parameters; u 0 , E 0 , 
7 , B, and T 0 . The last one is not very relevant; as long as it is chosen as a low temperature 
where the desorption is negligible the TPD spectra will be the same. We can remove two 
more parameters by scaling the time and the temperature. Instead of time we use i/ 0 f, and 
instead of temperature we use k B T/E 0 . The only remaining parameters are then 7 and 
k B B/ vqEq, the scaled version of B. The former is, of course, the more relevant one. We 
will also present results for desorption at constant temperature. These depend on just one 
parameter; k B T / E 07 . 

Fig. |T] shows the TPD spectrum for 7 = 0.05 and k B B/ vqEq = 10~ 16 . (This is roughly 
the TPD spectrum one would get with u 0 = 10 13 s -1 , E 0 /k B = 3000 K, E^/k B = 150 K, and 
B = 3K.s^, values very similar to the Xe/Pt(lll) system.0) Although 7 is rather small, 
we see a substantial shift of the peak maximum with decreasing initial coverage. Lower 
coverages may show larger desorption rates. This implies that the order of the desorption is 
negative, which has been observed in Xe/Pt(lll) experimental spectra before.0’0 Smaller 
values of the parameter 7 show smaller distances between the peak maxima, but not until 
7 is clearly smaller than 0.01 the curves of different initial coverages stop crossing. (We will 
derive a better estimate below.) Making 7 larger spreads out the peak as far as we could 
check. However, for 7 larger than about 0.1 we ran into convergence problems with the grid 
size and the diffusion rate. For really large values of 7 one expects island formation. Islands 
are one explanation for zero-order desorption ,0 so the trend towards more negative order 
should reverse for some 7 . 
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As the interactions are attractive between the adsorbates, there is a order-disorder phase 
transition. Below the phase-transition temperature islands just mentioned are formed. This 
phase-transition temperature is about half E^/ks■ The negative-order TPD spectra have 
peaks well above this temperature, however. This implies that a model, which assumes that 
the adsorbates are randomly distributed over the substrate, might explain the peak shifts. 
To check this it seems more convenient to look at desorption at a constant temperature, as 
this reduces the relevant parameters to just one; ksT/Eo'y. The results are shown in Fig. [2]. 
The model and the MC results agree very well except for the relatively low temperatures 
at intermediate coverages. At those temperatures there is still some short-range order. The 
adsorbates tend to have somewhat more neighbors than in a random distribution so that 
the desorption is smaller than in the model. The model and the MC simulation show that 
for desorption at low constant temperature starting with a high coverage the desorption rate 
increases, reaches a maximum, and then decreases. The initial increase implies a negative 
order. At higher temperatures the order is positive for all coverages, and so we have another 
transition temperature. This temperature can be defined as the temperature where the 
order is zero at 6 = 1. 

We can derive this temperature from the model, or more generally, by assuming that 
for 0 ps 1 there are only adsorbates which still have all their neighbors or have only lost 
one. With N a the number of adsorbates in configuration a the master equation gives us the 
general result 


d{N)_ 

dt 


d_ 

dt 


Y, P <* N « 


Y, w ^Pp{N a 

a/3 


AW 


-E p nE w<n ' N li n> . 

0 n =0 


( 5 ) 


where is the number of adsorbates in configuration (3 with n neighbors, and is the 
transition probability per unit time for such adsorbates. Dividing by the number of sites 
this gives us for the desorption rate r at 9 ps 1 the expression 
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( 6 ) 


dO 

dt 

= [! — (! + Z)( 1 - 9)}W {Z) + Z(1 - 9)W (Z ~ 1) , 


where Z is the number of neighboring sites. The coefficients of W^ are the fractions of 
all adsorbates with n neighbors. From this we can fold when the negative-order desorption 
disappears by setting dr/dt = — d 2 9/dt 2 = 0 at 6 = 1, which gives us 


W^ _ Z 
WV-V = Z + l' 


( 7 ) 


As W^ oc exp(— 7iE 0 r y/k B T), this leads to 


k B T = 


Eoj 


( 8 ) 


ln(Z + 1) — In Z 

for the transition temperature. For Z — 4 we have k B T fa 4.48 E^, which is indeed far 
above the order-disorder phase-transition temperature. It’s also instructive to compute the 
order x at 9 = 1. We find that it equals for desorption at a constant temperature 


x = 


9 fdr' 
r\d9 


= 9 


(P9_ 

dt 2 


= 1 


Y—V 

[Tt) 

0 E^/k B T 


( 9 ) 


For the values used in Fig. [l] this becomes x = —9.96 at the peak maximum temperature of 
k B T/E 0 = 0.0379. Whether the negative order really shows up in a TPD spectrum depends 
on the positions of the peaks. For normal heating rates the peak maximum is at about 
k B T / Eq fa 1/30. This value should be lower than the temperature at which the negative 
order disappears. Hence for Z = 4 we fold negative orders in TPD spectra if 7 > 0.00744. 

We would finally like to remark on the relevance of the results just presented to real 
systems. As has already been mentioned negative orders in TPD have been seen for 
Xe/Pt(lll).IHHI3 As the effect seems to be large, it should be easy to find experimentally. 
However, this depends crucially on the dependence of the activation energy for the desorp¬ 
tion on the number of neighbors. Above we have assumed that the activation energy was 
additive. If we assume a dependence E 0 ( 1 + y^A^) then we fold from Eq. (|7|) 
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( 10 ) 


k B T = £q 7 ' 


Vz-Vz^i 


ln(Z + 1) — In Z 

for the transition temperature. For Z — 4 this gives us k B T = 1.20 Eqj, and with the peak 
maximum at k B T/E$ « 1/30 we must have 7 > 0.0278 in order to get negative orders in 
TPD spectra. We see that the effects are much weaker than for additive interactions. As 
attractive interactions are primarily found for physisorption, the additive activation energy 
may seem to be a good approximation. However, substrate-mediated effects may give large 
non-additive contributions .0 

Non-additivity need not be the only factor obscuring negative orders as we can demon¬ 
strate for Xe/Pt(lll). Activation energies for this system as a function of the number of 
neighbors have been determined using molecular dynamics (MD)Jil They are well approx¬ 
imated by E 0 /k B = 2851 K and 7 = 0.031. As Z = 6 , negative-order desorption should 
disappear only at 575 K, whereas peak maxima are found somewhat below 110 K. Large 
peak shifts as a function of initial coverage are expected. Indeed, if we use the same acti¬ 
vation energies for the diffusion, we find spectra very similar to Fig. [IJ However, such large 
shifts have neither been found experiment.ally, 0 nor with MD.0 The reason for this is that 
the lateral interactions are larger than the activation energies suggest. A value of 7 = 0.049 
is more appropriate for diffusion. This is found using the effective Xe-Xe potential of the 
MD study. If we use this value for the diffusion we obtain the spectra in Fig. [3]. 

We see that the shifts are small. There is a reasonable agreement between Fig. [| and 
the zero-order-like spectra of the MD study. Figure [| still shows slight negative order, and 
there are no extra peaks or shoulders. We think that these latter features in the MD study 
are probably artifacts caused by fits and extrapolations. The same holds for the oscillations 
in the order and the prefactor (Fig. 7 of the MD study0), which are particularly sensitive 
to numerical errors. 

It might seem contradictory that attractive lateral interactions, that in our model with 
randomly distributed adsorbates cause large peak shifts, reduce these shifts in Xe/Pt(lll). 
The reason for this is that in the model the lateral interactions only effect the desorption, 
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whereas the effect in Xe/Pt(lll) is dominantly on the structure of the adlayer. For the 
model the probability that an adsorbate has n neighbors equals 9 n ( 1 — 9) z ~ n (note that 
this does not depend explicitly on the temperature). For Xe/Pt(lll) the fraction of Xe 
atoms with many neighbors is larger than according to this expression, as desorption takes 
place not far above the order-disorder transition. The probability of a Xe atom having n 
neighbors and the short-range order are similar to what was found in the MD study. Figure |2| 
already showed the consequence. The curves of the model are always above the results of the 
simulations. At low temperature this means, that starting at high coverage, the desorption 
rate in the model increases more than in the simulations, which implies a more negative 
order, and hence larger peak shifts. 

IV. CONCLUSIONS. 

TPD spectra of systems with attractive lateral interactions may show large shifts of the 
peak maxima as a function of initial coverages. These shifts indicate negative-order desorp¬ 
tion. A simple model with random occupation of sites can explain them qualitatively. The 
same model can also be used to derive an exact analytical expression for a transition tem¬ 
perature above which negative-order desorption is not possible. This transition temperature 
may be higher than the order-disorder phase-transition temperature by about an order of 
magnitude. Non-additivity and the structure of the adlayer may decrease the peak shifts. 
Therefore more experiments are necessary to see if the negative-order desorption that occurs 
in the model system described above also occurs in real systems. 
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FIGURES 

FIG. 1. TPD spectra showing the desorption rate — dQ/d(vot ) as a function of the temperature 
ksT/Eo for 7 = 0.05 and heating rate ksB/voEo = 10 -16 . The dots are results directly from the 
MC simulations, and the lines connect MC results with the same initial coverage. The five curves 
are from left to right for initial coverage 0.262, 0.463, 0.644, 0.822, and 1.000. 

FIG. 2. Scaled desorption rate — exp^o/fogT) dQ/d(yot) as a function of coverage 6. The 
symbols are results directly from the MC simulations, and the lines are from the model with 
random distribution of adsorbates. Results are shown for five temperatures from top to bottom 
ksT/Eo'y = 6 , 4.48, 3, 1.6, and 0.8. The value of ksT/E 07 = 4.48 corresponds to the transition 
temperature above which negative-order desorption is not possible. 

FIG. 3. TPD spectra showing the desorption rate —dO/dt (in number of Xe atoms per second 
per site) as a function of the temperature (in K) for heating rate B = 1 K-s _1 . The dots are results 
directly from the MC simulations, and the lines connect MC results with the same initial coverage. 
The five curves are from left to right for initial coverage 0.598, 0.698, 0.806, 0.906, and 1.000. 
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